Nonperturbative Dynamical Theory and A Scheme for Nonequihbrium Transport 



o 
o 

(N 
< 

o 



I 



I 

o 
o 



> 

vn ■ 

^ : 

o 
i> 
o 



I 

o 
o 



X 



Jongbae Hong and Wonmyung Woo 
Department of Physics and Astronomy & Center for Theoretical Physics, 
Seoul National University, Seoul 151-747, Korea 
(Dated: February 6, 2008) 

We develop a nonperturbative dynamical theory (NDT) that is useful for treating nonequihbrium 
transport in a system with strong correlation. We apply our NDT to the single-impurity Anderson 
model in equilibrium to check its reliability by comparing with the results of numerical renormal- 
ization group method (NRG). We finally suggest a self-consistent loop to calculate the current in a 
lead-dot-lead system with Kondo coupling. 



One of the most problematic subjects in the field of 
theoretical condensed matter physics is the treatment of 
nonequilibrium transport [ij in a system with strong cor- 
relation. There is no established theory for treating the 
interacting system when it approaches nonequilibrium 
state, and a strongly correlated system requires nonper- 
turbative treatment. Moreover, methods based on the 
static nonperturbative theory, such as NRG[1, 0], are 
inappropriate for studying nonequilibrium transport, be- 
cause current is a dynamical quantity. Therefore, a possi- 
ble method for the treatment of nonequilibrium transport 
in a system with strong correlation is by employing the 
NDT, which has not been studied successfully thus far. 
Here, we report an NDT that is developed by fully utiliz- 
ing the dynamical nature of the Heisenberg picture. We 
first apply the new NDT to the single-impurity Anderson 
model in equilibrium; further, we extend the application 
to the nonequilibrium transport of a lead-dot-lead sys- 
tem. 

Transport phenomena in mesoscopic systems have re- 
cently attracted considerable interest in connection with 
nonequilibrium transport in a strongly correlated system, 
which is one of challenging and debated subjects in recent 
theoretical condensed matter physics. An example of a 
typical quantum system with strong correlation under 
nonequilibrium conditions is a quantum dot with metal- 
lic leads under bias. Theoretical studies on quantum dot 
have usually reported the conductance properties that 
can be obtained by using the NRG, since conductance is 
an equilibrium property of the system. Even though the 
NRG has been successful in providing low-energy eigen- 
values and eigenstates rigorously for a part of strongly 
correlated systems, it cannot be a resolver of the nonequi- 
librium transport problem. 

A good formalism for nonequilibrium transport has 
been established by Meir and Wingreen^] in terms of 
nonequilibrium Green's functions. Let us consider the 
motion of a spin-up electron in a quantum dot under bias. 
Since the movement of a spin-up electron will be affected 
by the movement of the spin-down electron due to strong 
correlation at the dot, the retarded on-site Green's func- 
tion of a spin-up electron must contain information on the 
back and forth movement of the spin-down electron. A 
static theory such as NRG is incapable of yielding infor- 
mation on the back and forth movement of an electron 
with a particular spin. A dynamical theory, however. 



can incorporate this information in Green's function. In 
this work, we present an NDT giving the nonequilibrium 
Green's function that can provide a spectral density for 
the Anderson model when the system is in equilibrium, 
which is comparable to the one obtained by the NRG, 
and a scheme providing the current-voltage characteris- 
tics under nonequilibrium conditions. 

The Heisenberg picture stresses the dynamics of op- 
erators compared with other two pictures. However, 
the advantage of its dynamical nature has not received 
sufficient appreciation. Here, we focus on the dynam- 
ical nature of the Heisenberg picture to develop the 
NDT. The formal solution of the Heisenberg equation 
for a fermion annihilation operator Cd], i.e., Cd|(i) = 
Crft + [H, Cd^]{it) + [H, [H, Cd|]](it)V2 -I- • • • , where H is 
the Hamiltonian of the system, describes the time evolu- 
tion of Cd] ■ The operators in each term of the expansion 
represent linearly independent ways of annihilation of a 
spin-up electron at site d at time t. These operators pro- 
vide linearly independent bases spanning the operator 
or Liouville space. We call these the dynamical bases. 
The inner product in the Liouville space is defined as 
{A\B) = {{A,B'''}), where A and B are elements of the 
Liouville space, the curly brackets denote the anticommu- 
tator, and the angular brackets indicates an equilibrium 
or nonequilibrium average, depending on the situation. 
Constructing appropriate dynamical bases is the essence 
of the NDT. 

We first develop our NDT for the single-impurity 
Anderson model represented by the Hamiltonian 

^kd^lcr^da) + UndiUdi, whcrc Udi = 4|Cdx, and the sub- 
scripts fj, fc, and d denote the spin index, quantum state 
of the metallic reservoir, and the position of impurity, 
respectively, and then compare the result for the spec- 
tral density with that of the NRG. Finally, we apply our 
NDT to study the nonequilibrium transport of the lead- 
dot-lead system with on-site Goulomb interaction. 

The basis operators created by the commutators 
with the Hamiltonian are composed of Ck] where k = 
d, 1, 2, • • • , cx) and the operators combined with Ck] and 
other operators such as n^x, jd], and j'^^, where j^^ = 

i[H,ndi] = iiJ2k VfcdcIiCdi-Efe 14d4iCfei)' i-^-, the cur- 
rent operator, and j+ = i[ndi,jdi\ = CLk^kdcl^Cdi + 
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FIG. 1: Pictorial descriptions of the roles of basis operators 
UdiCk^ (a) and j^|Cd| for a particular k-state (b). Sum over all 
k in (b) becomes j^f^Cdf. The dashed arrows do not indicate 
hybridization. They just indicate creation and annihilation. 



J2k^kd^di^ki)- We keep the meaningful operators that 
describe the virtual exchange between spin-up and spin- 
down electron and take the mean-field approximation for 
other part of the primitive basis operator. After this ma- 
nipulation, appropriate linearly independent dynamical 
bases spanning a reduced Liouville space of the operator 
Cdi can be constructed. The bases except Cd-\ itself must 
be orthogonal to c^i in order to yield a correct projec- 
tion ({c|j|, C(i|(i)}) that gives the on-site retarded Green's 
function. 

The dynamical bases of the reduced Liouville space of 
Cd] are composed of five parts. The first two are (i) a 
set of bases Sk defined by Sk = {ck^\k = l,2,---,cx)} 
for describing the annihilation at the impurity site af- 
ter some hoppings in the metal and (ii) a set of bases 
Sn defined by Sn ~ {SridiCk-ilk = 1,2,- ••,00} where 
Sudi — ndi — (ndi) for describing the number fluctuation 
of the spin-down electron at the impurity site during the 
annihilation process of (i). Lastly, we consider the bases 
coupled to Crf|, i.e., (iii) a set of bases Sd defined by 
Sd = {cdh^jdi'^dh^Jdi'^d]} for describing the annihila- 
tion of a spin-up electron at the impurity site without 
any coupled processes and with the processes coupled to 
the back and forth fluctuation of the spin-down electron 
at the impurity site. The latter describes the hybridiza- 
tion in the Kondo process. 

The bases of the first subset Sk and c^f constitute a 
complete set when the system is noninteracting. Since 
Sn'^^ is composed of a constant and Sndi, the first two 
subsets, Sk and Sn, along with c^f, contain all the dy- 
namical processes described by number fluctuations of a 
spin-down electron during the annihilation process of a 
spin-up electron at site d. The operators in the third 
subset, Sj^^Cdi, contribute to describing the Kondo ef- 
fect. As an example, the pictorial descriptions of the 
basis operators SudiCki and Sjf^Cdi are shown in Fig. 1. 
We do not consider the other complicated processes that 
can occur during the annihilation process of a spin-up 



electron at time t since the processes considered above 
may be sufficient for describing the Kondo phenomena in 
the single-impurity Anderson model. In addition, the op- 
erator Cdi is orthogonal to all other bases of the reduced 
Liouville space. 

We are now in a position to obtain the spectral den- 
sity of a spin-up electron at site d by NDT using the dy- 
namical bases introduced above. The resolvent Green's 
function operator in the Heisenberg picture is written 



as 



(w ± iry — L) , where L is the Liouville op- 



erator defined by hA ~ HA — AH, 77 is a positive in- 
finitesimal, and the superscripts ± denote retarded (+) 
and advanced (— ), respectively. This expression becomes 
— {uj ± irj — H)~^ in the Schrodinger picture. 

Since the retarded and advanced Green's functions are 
given by Gfj{oj) = (ei|(a; ± ir] - L)"i|ej), where \ej) 
is one of basis operators spanning the Liouville space, 
they are given by the inverse of the matrix M defined 
by = (eijzl -I- iL|ej), where z = —iuj ± 77, i.e., 

iGfjiuj) = (e,|M"i|ej) = (adj M)y[det M]-\ where 
(adj M)y denotes the cofactor of the jz-element in the 
determinant of M^Sj]. This expression has also been re- 
ported in literature [6|. 

In constructing matrix M, we use normalized bases in 
order to make it quasi-symmetrical. If we arrange Sk, Sn, 
and Sd in a regular sequence to construct the matrix M 
for the single-impurity Anderson model, it is represented 
by a matrix of four blocks. 



Ms 



lAM 



Mr Mdr 
-M*dr Md 



where the subindex r represents metallic reservoir. The 
blocks Mr, Mdr, and M^ are 00 x 00, 3 x 00, and 3 x 
3 matrices, respectively. Since the inner product has a 
relation {ihA\B) = —{ihB\A)*, only the real parts of 
the matrix elements have different signs for their counter- 
parts. 

The block M^ is represented by 



Mr 



Mil 
Mil 



where Mn is diagonal and its elements are z + ie^, k = 
1, 2, • • • , CJO. The 3 X 00 block Mdr, on the other hand, 
has the form 



M, 



Ckd 








where the column Ckd has iVkd as its elements, while 
Cnj- and C„j+ have idVkd and £,'^Vkd, respectively. 
The blocks except M^ contribute to the self-energy. 
The 11-element of the block M^ is 2: -f- i{ed + {ndi)U), 
and the other diagonal elements are given by 2: -f icd + 
\U {ndi5j^^) / {dj"^^)]. These are equal to the 11-element 
of M-d under the decoupling approximation. The off- 
diagonal elements of M^, on the other hand, are given 
by = -Mf = (C//2)C, = -Mf = (£7/2)^+, 
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and Mf 



-Mf 



7, where 2^^ = - 2(ndi))(jl) + 
(zK,jJJ(l~2n,T))]/[((5jl,)^)i/2((<5n,j2)V2] and 7 = 

-2*Efc ^*.07ijicfcT4T)/[(('5M)'>(('5ji)')]'/'- These 
factors will be determined when we calculate the spec- 
tral densities. 

In order to handle the infinite dimensional matrix, 
matrix reduction by Lowdin's partitioning technique^ 
is performed by solving the eigenvalue equation for the 
original matrix Msiam, such as MsiamC = 0, where C 
and are infinite dimensional column vectors. The col- 
umn vector C is partitioned into two parts, i.e., infinite 
dimensional and three dimensional Cd, symbolizing 
the reservoir and impurity parts, respectively. Then, the 
equation for is obtained as {'M.d—Mrd^d^^dr)Cd = 
MrfCc; = 0. The reduced 3x3 matrix contains 
the information on the many-body dynamics of a spin- 
up electron starting from the impurity site at t = and 
ending at it at time t. 

It is impossible to calculate the inverse of a general 
(XI X cx) matrix; however, this is not the case for the matrix 
Mr, which is block diagonal. The inverse of can 
be obtained in a straightforward manner The second 
term appears as additional self-energy terms in after 
reduction. One can easily imagine that = ^j" in the 

Kondo regime at half-filling, and the final form of for 
a symmetric Anderson model in which = —U/2 and 
(rid) = 1/2 is given by 



Md = 



( —iuj + 



7 + ^if^o 



UCdl2 \ 



iuj + i^^ ^So / 



where Eo = Ek |VfcdP/(w - ek + iv) = H^) - «^(^) 
is the self-energy of the Anderson model with U = 0. 
We will neglect the k dependence of Vkd in this work. 
The real and imaginary parts are respectively A(w) = 
(P/tt) A{€)de/{uj — e), where P represents the prin- 
cipal integration and A(a;) = TT^^.\Vkd\'^^{^^ — Ck) = 
A[l — (w/Z?)^]^/^ for the semielliptical band, where A = 
2\V\'^/D and D is half of the band width. Now, one 
can obtain the retarded Green's function that is given 
by iG'^j^{uj) = (adj Md)ii[det M^]"^ by calculating the 
inverse of matrix M^. 

We find that 7 determines the width of the Kondo 
peak, while governs the spacing of the side- 
peaks of the spectral density. Therefore, from the 
analysis at atomic limit, one can find that = 
I/V2. On the other hand, our retarded Green's func- 
tion gives the real and imaginary parts of the self- 
energy as ReE^jjiuj) = -{U'^ /4j'^)uj + 0{uj^) and 

(l/A)ImI]+;(w) = {U^/2-f'^)L0^+Oiuj'^). The subscript 
indicates the interaction part of the self-energy. The real 
part of the self-energy provides the wavefunction renor- 
malization as Z =1/[1 + {U^ /i'j'^)]. 



FIG. 2; Spectral densities of the symmetric Anderson model 
with various Coulomb repulsions. 



Unfortunately, it is difhcult to obtain 7 by direct cal- 
culation. Therefore, we determine rigorous 7 by using 
the result of Bethe ansatz, which is given by Z^^ — 
{4:/T:)y/U/rcxp[-TTU/4T+TTr/4U] for the symmetric An- 
derson model[3]. We plot the spectral density Pd^{i^) = 
— (l/7r)ImG'J^(tj) in Fig. 2. Our result naturally com- 
patible with that of NRG at least in the Kondo regime. 

The above result can be extended to nonequilibrium 
transport in a lead-dot-lead system with on-site Coulomb 
interaction at the dot in a straightforward manner be- 
cause it is simply a single-impurity Anderson model with 
two separate metallic reservoirs. This only necessitates 
the simple requirement of an additional number of bases 
for describing the left and right leads and the movements 
from or to both leads, such as 5*^, Sf^ for the left lead, 
5^ = {Sjdtcd]'^Sj+l'cd^,Cd^,6j-^''cd^,Sj+^''cd]} for the 
dot, and Sj}, for the right lead. The superscripts L 
and R denote the left and right leads, respectively. 

If we arrange 5*^, S*^, S^"^^, S^, in a regular se- 
quence to construct the matrix M for a lead-dot-lead 
system, it will be represented by a matrix of nine blocks. 














Mid 




MLd 




M.Rd 






K 




Mrr 



where blocks M^, M^^ and M^ij, and M^^ and M^^ are 
5 X 5, 5 X 00, and 00 x 5 matrices, respectively. Blocks 
and are 00 x 00 matrices that are constructed 

by the sets of bases describing left and right leads, re- 
spectively. Since no direct coupling exists between the 
left and right leads, zero matrices are present at two of 
the corners. The structure of each block is similar to that 
of the corresponding block of the matrix Msiam. 

The matrix reduction produces a 5 x 5 matrix for M^j, 
whose inverse yields the on-site retarded Green's func- 
tion G'^^{uu) of the lead-dot- lead system in the frame- 
work of the NDT, which is given by a function of (uda), 
Uda"'^)^ Uda '^)- Siucc the Current is expressed by 
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J - 



- J« = {J!;-J^)I2, where = {e/h)U^J^), 
the current and other variables are expressed by the lesser 
and retarded Green's functions such that [J, [j] 



07.^) = J -{^^'A^)G^M 

- 2/i(c.)r^(c.)ImG+,(a;)} 



(1) 



Or)= / ^[/LHr^(c.')ReG+,(c.)-*G<,,(a;)A(, 



(2) 



(3) 

where the effective Fermi distribution /{uj) is given hy^ 



In order to construct a self-consistent loop for cal- 
culating the current-voltage characteristics, we use the 
Keldysh equation[i| 



G<icj)^G:iu;)J:<{cj)G-iu;). 



(5) 



Then, one can construct a self-consistent loop from Eqs. 
(l)-(5) as 

This iterative method for calculating noncquilibrium 
quantities of the strongly correlated system is the sec- 
ond major result of this work. The specific results will 
be reported in a separate work. 
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